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ABSTRACT 


A Type Ia supernova (SN) at z = 1.78 was discovered in James Webb Space Telescope Near Infrared 
Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; z = 0.35). The SN is situated 1.5- 
2kpc from its host galaxy Arc 2 and appears in three different locations as a result of gravitational 
lensing by G165. These data can yield a value for Hubble's constant using time delays from this 
multiply-imaged SN Ia that we call “SN HOpe.” Over the entire field we identified 21 image multi- 
plicities, confirmed five of them using Near-Infrared Spectrograph (NIRspec), and constructed a new 
lens model that gives a total mass within 600 kpc of (2.6 + 0.3) x 101^ Mọ. The photometry uncov- 
ered a galaxy overdensity at Arc 2's redshift. NIRSpec confirmed six member galaxies, four of which 
surround Arc 2 with relative velocity X900 km s-! and projected physical extent <33 kpc. Arc 2 
dominates the stellar mass ((5.0 + 0.1) x 10!! Mo), which is a factor of ten higher than other members 
of this compact galaxy group. These other group members have specific star formation rates (sSFR) 
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of 2-260 Gyr-! derived from the Ha-line flux corrected for stellar absorption, dust extinction, and slit 
losses. Another group centered on the dusty star forming galaxy Arc 1 is at z = 2.24. The total SFR 


for the Arc 1 group (2400 Mo yr^!) translates to a supernova rate of ~1 SNe yr 


-l, suggesting that 


regular monitoring of this cluster may yield additional SNe. 


Keywords: large-scale structure of universe - gravitational lensing: strong — galaxies: fundamental 
parameters — galaxies: clusters: general — galaxies: high-redshift 


1. INTRODUCTION 


The galaxy cluster PLCK G165.7+67.9 (G165; 
z = 0.35) first garnered attention by gravitationally- 
amplifying a single, apparently-bright galaxy in the 
background called *Arc 1." Arc 1 is boosted to an 
observed sub-mm flux density of S35o,;,, > 700 mJy 
(Canameras et al. 2015; Harrington et al. 2016), making 
it detectable by the Planck and Herschel Space Observa- 
tory missions (Planck Collaboration et al. 2016, 2020a). 
In addition to its lensing properties, G165 also has some 
other unusual physical properties. It is a double-cluster 
with two prominent cluster cores, its X-ray luminos- 
ity is lower than galaxy clusters of a similar redshift, 
and its mass, based on the Planck SZ Compton-Y map, 
yields only an upper limit. These properties arise in 
part because the cluster has a relatively-small mass of 
2—3x1014 Mo (Frye et al. 2019; Pascale et al. 2022), and 
in part because the source responsible for the high rest- 
frame far-infrared flux density is not the entire galaxy 
cluster but a only portion of a single lensed source that 
is a dusty and high star forming galaxy (DSFG, Casey 
et al. 2014) that we call herein *Arc 1 b/c." 

Strong lensing by G165 renders Arc 1 into two images 
merging with the critical curve, Arc 1b/c and a sepa- 
rate counterimage that we call herein “Arc la,” which 
together form an “image system". The star formation 
rate (SFR) for Arc 1b/c uncorrected for lensing mag- 
nification was measured to be uS FR ~ 12000 — 24000 
Mo yr} that was based on integrating the spectral en- 
ergy distribution (SED) fit over a wavelength range of 8- 
1000 um (Harrington et al. 2016). A requisite byproduct 
of this high star formation is UV radiation by massive 
stars. Interestingly, Arc 1b/c is detected in the rest- 
frame UV (observed frame g-band) (Frye et al. 2019). 
'This is impactful because the vast majority of the light 
is expected to be obscured by dust (Casey et al. 2017; 
Whitaker et al. 2017). The UV light may be made more 
detectable by the high lensing magnification factor of 
Z, 30 that spreads the light over more pixels, thereby 
potentially enabling higher signal-to-noise detections of 
the escaping UV radiation field. The physical conditions 
are less well constrained owing to the lack of a rest-frame 
optical spectrum with which to obtain the galaxy clas- 


sification (i.e., Mingozzi et al. 2023), estimate the dust 
extinction, investigate the star-forming properties via 
the Balmer emission lines (i.e., Kennicutt 1998; Shap- 
ley et al. 2022), and the measure the gas-phase oxygen 
abundance (i.e., Curti et al. 2020; Li et al. 2023). 
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Figure 1. JWST/NIRCam coverage of the G165 field. The 
background is the r-band negative image from Megaprime. 
Color images show the combined NIRCam data. The pink 
long-dashed rectangle outlines Epoch 1, and Epochs 2 and 3 
are squares that mostly overlap but have slightly different 
rotation angles. The blue square outlines the field of view 
of previous HST WFC3-IR imaging, which usefully covers a 
portion of the gap between the two NIRCam modules. The 
green square outlines the field of view adopted to construct 
the lens model. 


Another prominent infrared-bright galaxy is Arc 2 
(mg = 20.5 AB mag), which has a robust photometric 
redshift based on a joint fit of the imaging from seven 
ground- and space-based facilites of z = 2.30 + 0.32 
(Pascale et al. 2022). This redshift has recently been 
confirmed spectroscopically from the ground at zsp = 
1.783 + 0.002 (Polletta et al. 2023). Somewhat surpris- 
ingly, in Pascale et al. (2022) the photometric redshifts 
for several arc systems hovered around the same red- 
shift, zph ~ 2. It is tempting to conclude that there is a 
galaxy overdensity at this redshift, but existing imaging 
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was inhomogeneous, with limiting magnitudes ranging 
from 23.3-28.9 AB and PSF point-source full-width-half- 
maximum values ranging from 0.13-2.02", complicating 
their confirmation as z ~ 2 sources by their photometric 
redshifts alone. 

HST WFC3-IR imaging enabled the identification of 
large sets of image systems and the construction of 
lens models (Cafiameras et al. 2018; Frye et al. 2019). 
These models confirmed both the northeastern (NE) and 
southwestern (SW) cores of this binary cluster, and the 
high, cluster-scale mass. This earlier work was built and 
expanded upon by the addition of robust photometric 
redshift estimates across ground and space-based data 
sets for three of the 11 image systems and several cluster 
members (Pascale et al. 2022). Yet, these models were 
anchored on the spectroscopic redshift for only one im- 
age system (Arc 1), thereby limiting the accuracy of the 
resulting lens model and its ability to recover the lensed 
image positions (Johnson & Sharon 2016). 

On account of its physical properties, prominently- 
placed infrared galaxies, and ample strong lensing ev- 
idence, G165 was observed by the James Webb Space 
Telescope (JWST) as a part of the Prime Extragalactic 
Areas for Reionization and Lensing Science (PEARLS) 
program (Windhorst et al. 2023). The PEARLS Near 
Infrared Camera (NIRCam) imaging, referred to as 
“Epoch 1," uncovered three point sources which were 
bright (mrisow = 23.91+0.01 AB mag) yet not present 
in the previous HST imaging. This triply-imaged tran- 
sient was estimated to be a Type Ia supernova (SN) with 
at a > 90% probability based on the three points popu- 
lating the light curve obtained in Epoch 1. A JWST 
disruptive Director's Discretionary Time (DDT) pro- 
gram was executed that contained two additional imag- 
ing epochs which sum up to a total of 9 samplings of the 
supernova SN light curve (three per image per observ- 
ing epoch). Epoch 2 imaging was scheduled three weeks 
after Epoch 1, so that a set of photometric points could 
be taken just prior to the second peak of the SN Ia light 
curve. Epoch 3 imaging was scheduled two weeks after 
Epoch 2 to bracket the second SN peak on the fading 
side. NIRSpec spectroscopy was obtained for all three 
SN appearances during Epoch 2 when they were rela- 
tively bright, and included the SN host galaxy Arc 2, 
and other high-priority arcs (PID: 4446, PI: Frye). 

We present herein an overview of the initial science re- 
sults in the combined JWST PEARLS and DDT imag- 
ing and spectroscopic observations in the G165 cluster 
field and of the discovery of the SN, which is dubbed 
“SN HOpe” for its potential to measure the time delays 
between the images, and from it, a value for Hubble's 
constant (i.e., Kelly et al. 2023a,b; Liu et al. 2023, and 


references therein). This study is the first in a series of 
papers whose objective is to investigate SN HOpe, the 
cluster, and the lensed sources. 

'This paper is organized as follows. 82 introduces the 
JWST data sets and the ancillary data. $83 describes 
the NIRCam photometry, the photometric discovery of 
SN HOpe, and the estimation of photometric redshifts. 
The construction of the lens model appears in §4. The 
NIRSpec spectroscopic analysis follows in 85 with a fo- 
cus on the z z 2 galaxy groups. The physical properties 
of these galaxy groups are investigated in 86. 87 sum- 
marizes the results. This paper uses the AB magnitude 
system throughout, and we assume a flat ACDM cos- 
mology with Ho = 67 km s^! Mpc™t, Om,o = 0.32, and 
Qao = 0.68 (Planck Collaboration et al. 2020b). 


2. OBSERVATIONS AND REDUCTIONS 
2.1. NIRCam 


'The Epoch 1 NIRCam observations were obtained as 
part of the PEARLS program (PID 1176, PI: Wind- 
horst). The observing window was selected by STScI 
in order to reduce stray light expected from a nearby 
bright star. Exposures were taken in four filters in the 
short wavelength (SW) module and four in the long 
wavelength (LW) module as shown in Table 1. Both 
NIRCam modules collected data. Epochs 2 and 3 of 
NIRCam imaging were acquired as part of a JWST dis- 
ruptive DDT program (PID 4446, PI: Frye) to follow 
the supernova's light curve in each of its three images. 
In this follow-up program exposures were taken in six 
filters using only Module B of NIRCam, also recorded 
in Table 1. The NIRCam observations covered the cen- 
tral region of the cluster including both the NE and 
SW cluster components, the three images of the SN, the 
DSFG, all of the image systems, and the prominent gi- 
ant arcs. Figure 1 depicts the field coverage relative to 
the Megaprime r-band image. 

The images were reduced by our team as described 
by Windhorst et al. (2023). Briefly, the data were re- 
trieved from the Mikulski Archive for Space Telescopes 
(MAST), and the latest photometric calibration files 
were used (pmap.1100). All images were reduced using 
version 1.11.2 of the STScI JWST Pipeline (Bushouse 
et al. 2022) with an additional correction for 1/ f noise 
by applying the prescription of C. Willott!. The Pro- 
Found code was run, which makes a second round of 
corrections of other image artifacts in the relevant rows 
and columns. This step additionally flattens the back- 
ground and corrects for detector-level offsets, “wisps,” 


1 https://github.com/chriswillott /jwst.git 
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Figure 2. NIRCam color image of the central region of G165. G165 is a double cluster with prominent NE and SW components. 
Colors follow the prescription in Trilogy (Coe et al. 2012) with red showing F444W and F356W, green showing F277W and 
F200W, and blue showing F150W and F090W. The 21 image systems used in our lens model are marked. They include the 
DSFG as Arcs la and 1b/c. The triply-imaged SN is labeled as “SN 2a/2b/2c." The orientation and image scale are provided 


for reference. 


and “snowballs” (Robotham et al. 2017, 2018). Since 
the Windhorst et al. (2023) publication, improvements 
in the data reduction techniques have been made by 
Robotham et al. (2023) regarding the removal of im- 
age wisps by using the wisp-free LW images as priors to 
identify the outer contours of the real detected objects. 
Those real objects were subsequently removed from the 
SW images to get a pure wisp image which was then 
fully subtracted. This process yields an image noise in 
the final mosaic that is almost the same in the wisp- 
removed area as in the surrounding wisp-free areas. 
After each frame was calibrated, the frames were 
aligned onto a common astrometric reference frame and 
drizzled into mosaics with pixel scale 20 milli-arcseconds 
(mas). The process was similar to that first described by 
Koekemoer et al. (2011) but updated to use the JWST 
pipeline”. Mosaics were produced for each filter in each 
separate epoch. For the six filters in common, all epochs 
were combined into a grand mosaic for each filter. AII 


? https:/ /github.com/spacetelescope/jwst 


mosaics were aligned onto the same pixel grid based on 
deep, ground-based Megaprime images with good see- 
ing, which had already been aligned directly onto Gaia 
DR3. The NIRCam data were aligned onto this grid 
with net offsets below 2-3 mas and no significant large- 
scale distortions. Figure 2 shows the central region of 
G165 in the main NIRCam mosaic. 


2.2. NIRSpec 
2.2.1. Data Reduction 


NIRSpec medium-resolution Micro-Shutter Array 
(MSA) spectroscopy of the G165 field was obtained 
on 2023 Apr 22. The MSA mask was populated 
with the positions of the three SN appearances (SN 
2a, 2b, and 2c) and two of the three images of the 
SN host galaxy (Arc 2a and 2c), followed by coun- 
terimages of three other image systems (Arcs 5a, 8c, 
and 9c). The remainder of the mask was filled with 
other lensed sources which summed to a total of 42 
lensed targets. The observations used the grating/filter 
combinations G140M/F100LP to cover spectral range 
0.97-1.84 um (rest-frame 0.35-0.66 um at z = 1.8) 
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and G235M/F170LP to cover 1.66-3.17 um (z = 1.8 
rest-frame 0.57-1.1 um), both with spectral resolution 
R = 1000. We also acquired a PRISM/CLEAR spec- 
trum covering 0.7-5.3 um (rest-frame 0.25-1.9 um) 
with R ~ 20-300 (50-14 A). All of the seven sup- 
plied guide stars were acquired, resulting in especially 
tight pointing residuals of 1-7 mas and successful point- 
ing even for targets near the edges of an MSA array. 
The science exposure times were 4420 s, 6696 s, and 
919 s for G140M/F100LP, G235M/F170LP, and the 
PRISM/CLEAR observations, respectively. A 3-point 
nod pattern was selected for each observation, and each 
MSA slit consisted of 3 microshutters giving slit height 
1.52". MSA slits are 0720 wide in the dispersion direc- 
tion, and the long dimension was oriented at position 
angle 276?. 

The Stage 1 calibrated data were retrieved from 
MAST, and reduced using the JWST NIRSpec pipeline, 
version 1.11.3.? Stage 2 and 3 reduction used the JWST 
pipeline with reference files “jwst_1100.pmap” for all lev- 
els of the data reduction, with an exception regarding 
the background subtraction, as is described below. Satu- 
rated pixels and other image artifacts were flagged in the 
2D spectra. The NIRSpec IRS2 detector readout mode 
was used, which largely reduced the 1/ f noise. The 2-D 
spectra were wavelength- and flux-calibrated based on 
the calibration reference data system (CRDS) context. 
Finally, individual calibrated 2D spectra exposures were 
coadded, and one-dimensional (1-D) spectra plus uncer- 
tainties were optimally extracted (Horne 1986). 

'The pipeline background subtraction performed well 
for single point sources and single small sources which 
were fully covered by the aperture. This is because 
the dithered exposures provided a good “best-fit” back- 
ground consisting of the intracluster light and/or other 
underlying extended sources, and/or detector noise lev- 
els. Hence, the resulting NIRSpec flux from the pipeline 
directly gave the flux for the point/small source. 

However, we lacked exposures of a separate back- 
ground target/field, which made it more of a challenge 
to estimate the background for sources extending across 
multiple microshutters. One example is the SN host 
galaxy Arc 2 and the SN, for which all three microshut- 
ters are occupied by sources. For this case the back- 
ground template formed in the NIRSpec pipeline comes 
from the flux through the source shutter in the dithered 
exposure. This “image from image" background might 
include some flux from the galaxy, leading to an oversub- 
traction of the background. A complementary problem 


3 https://jwst-pipeline.readthedocs.io/en/latest /index.html 


Table 1. JWST Epochs and NIRCam Ex- 
posure ‘Times 


Filter Epoch 1 Epoch 2 Epoch 3 
Mar 30 Apr 22 May 09 


F090W 2491 1246 1417 
F115W 2491 Zr ee 

F150W 1890 859 1246 
F200W 2104 1761 1761 
F277W 2104 1761 1761 
F356W 1890 859 1246 
F410M 2491 tee ee 

F444W 2491 1246 1417 


NoTE—Dates of each epoch are in 2023, and 
exposure times for each filter are in seconds. 


is the case for which neighboring microshutters are occu- 
pied by different sources. In the “MOS Optimal Spectral 
Extraction” tool from STScI (based on the method from 
Horne (1986)) a source kernel and a polynomial back- 
ground template are fit at the same time for one source 
within an MSA slit, based on a spatial window in the 2D 
spectrum chosen manually by the user, but the software 
does not support multiple source extraction. 


To alleviate some of these issues, a custom-built code 
was developed to perform the background subtraction. 
The code is different from the pipeline in that it builds 
a more locally-derived background template. For each 
pixel we evaluated the minimum flux of the set of five 
dithered pixels. Then for each pixel within each spatial 
column i the best value for the background was com- 
puted by the median value of this minimum flux within 
a running boxcar centered on column 7 that is 10 spa- 
tial columns wide. We found 10 columns to be a good 
compromise, for which a smaller median filter starts to 
encroach on the size of a typical cosmic ray mask, and 
a larger median filter smoothed out the background fea- 
tures in this wavelength-dependent operation. To cope 
with image crowding, the code has a multiple source 
extraction mode that fits multiple source kernels simul- 
taneously for each object along the MSA slit. 

Operationally, we ran NIRSpec Stage 2 with the back- 
ground subtraction task turned off, and then applied the 
custom-built code. The detailed content of this code, 
and its implementation for this data set, appear else- 
where (Chen et al. 2023). The results are recorded in 
Table 2. 


2.2.2. Spectroscopic Results 
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Figure 3. NIRCam color image centered on one of the 
SN host images 2c (Arc 2c), along with three other galaxies 
spectroscopically-confirmed to be at z = 1.78 in this study. 
Our lens model predicts for Arcs 2, 8, 9, and NS_337 to be sit- 
uated within 33 physical kpc in the source plane. The NIR- 
Spec spectra of these arcs appear in Figure 4. Colors follow 
the prescription in Trilogy (Coe et al. 2012) with red show- 
ing F444W and F356W, green showing F277W and F200W, 
and blue showing F150W and F090W. 


In total 47 sources were extracted from the NIRSpec 
dataset. We measured the redshifts from emission- and 
absorption-line features of each source, as available. The 
line centers were determined by fitting Gaussians to each 
spectroscopic line feature using SPECUTILS (Earl et al. 
2023). The 30 of 47 sources with secure measured red- 
shifts are listed in Table 2. (A redshift is considered se- 
cure if it has a high-significance detection of two or more 
spectral features and > 2c level in the continuum.) Of 
these, the highest-redshift galaxy is NS_274, a relatively 
rare example of a quiescent galaxy at z = 4.107, and 
the highest redshift multiple image is Arc 5.1a with a 
redshift measured off of Balmer lines from Ha through 
He detected in emission of z = 3.9530 + 0.0004. We 
note thaat the prefix “NS” stands for NIRSpec, and pre- 
cedes all of the NIRSpec-confirmed galaxy images in this 
study. 

The three SN HOpe spectra are of high quality. Most 
prominent is the requisite detection of the Si II \6355 
absorption feature blueshifted to ~6150 A, closely fol- 
lowed by the detection of the [Ca II] \AA8498,8542,8662 
IR triplet “CaT,” amongst other spectroscopic features. 
The spectrum of the SN, the SN classification as Type 
Ia, and the measurement of the spectroscopic time de- 


lay will appear in a different paper (Chen et al. 2023, in 
prep). The spectra of the SN host galaxy Arcs 2a and 2b 
had pre-existing redshifts, both based on the joint detec- 
tion of [O II] \3727 and the 4000 Å and Balmer breaks 
(Polletta et al. 2023). The new NIRSpec spectra give 
the first redshift for Arc 2c, whose value matches that 
of Arc 2a (Table 2). Nearly 20 spectroscopic features 
are detected as well as the 4000 À and Balmer breaks. 
Somewhat remarkably, six different lensed sources have 
the same redshift as the SN. Their images are shown in 
Figure 3 and their spectra are presented in Figure 4. The 
line fluxes and equivalent widths for several absorption 
lines in Arc 2 are presented elsewhere (Chen et al. 2023, 
in prep.). This redshift is interesting because it also co- 
incides with the strongest peak in the photometric red- 
shift distribution after the cluster redshift, as described 
in § 3.3. 

Arc 1b/c had a previously measured z,, = 2.2357 + 
0.0002 (Harrington et al. 2016). The NIRSpec spec- 
trum (NS_969) gives the first redshift of its counterim- 
age, Arc la, and the redshifts agree. A second spectrum 
(NS.46) in an adjacent MSA microshutter (0746 north- 
west) is redshifted by —1400 km s^! relative to Arc la. 
Both spectra, shown in Figure 5, have strong nebular 
emission lines and star-bursting attributes. Section 5.2 
describes their nebular and stellar properties. 

In some cases, most notably for Arcs 2a, 2c, and 
Arc NS.337, aperture losses are expected because the 
MSA slit coverage is smaller than the source size. To 
account for this shortfall, for each filter bandpass syn- 
thetic photometry was computed based on the NIRSpec 
PRISM spectra, which provided continuous coverage 
over the wavelength range of all eight NIRCam bands. 
We then compared our results to the NIRCam photom- 
etry integrated over the entire source. These correction 
factors are discussed in more detail in 85 and have been 
applied to the results in Table 4, where indicated. 


2.3. Other Optical and Near-infrared Observations 


LBT Large Binocular Camera gi-band imaging was 
acquired on 2018 January 20 (2018A; PI: Frye). Images 
reach 30 limiting magnitudes of 25.42 and 24.67 for g 
and i, respectively (Pascale et al. 2022). In the near- 
infrared, HST WFC3-IR exposures were taken on 2016 
May 30. Images at F110W and F160W reach 3o limit- 
ing magnitudes of 28.94 and 27.97 AB mag, respectively 
(Cy23, GO-14223, PI: Frye). Details of the observations, 
reduction, and analysis were given by Frye et al. (2019). 

At longer wavelengths, LBT LUCI-ARGOS imaging 
was obtained in K (2016B; PI: Frye). LUCI+- ARGOS 
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Table 2. NIRSpec Spectra 


ID R.A. Decl. ™MF200W ,obs Zsp 

NS.2 (SNa) 11:27:15.31  --42:28:41.02 5 b 

NS.3 (SNb) 11:27:15.60 — --42:28:33.73 a b 

NS_4 (SNc) 11:27:15.94 +42:28:28.90 2 b 
NS.6 (Arc 2c) 11:27:15.08 +42:28:28.72 20.26 1.7834 + 0.0005 
NS.7 (Arc 2a) 11:27:15.34 — --42:28:41.05 20.30 1.7833 + 0.0010° 
NS.19 (Arc 5a) — 11:27:13.20 +42:28:25.73 25.36 3.9530 + 0.0004 
NS.26 (Arc 8.2c)  11:27:15.89 +42:28:28.90 26.13 1.7839 + 0.0002 
NS_29 (Arc 9c) 11:27:16.11  +42:28:27.99 26.10 1.7816 + 0.0002 
NS.46 11:27:13.87 — --42:28:35.67 24.04 2.2401 + 0.0002 
NS_104 11:27:00.84 +42:27:03.65 24.83 3.1109 + 0.0004 
NS-_112 11:27:02.46  +42:27:10.22 21.59 0.6227 + 0.0003 
NS_123 11:27:04.64 +42:27:15.24 23.94 1.7874 + 0.0003 
NS_143 11:27:07.64 +42:27:26.22 23.28 1.6322 + 0.0002 
NS_171 11:27:05.82  +42:27:37.10 23.33 1.1787 + 0.0002 
NS_274 11:27:11.67 +42:28:10.70 23.44 4.1076 + 0.0023 
NS_285 11:27:06.82  +42:28:12.88 23.81 0.4466 + 0.0001 
NS_337 11:27:16.28 +42:28:23.59 21.77 1.7810 + 0.0009 
NS_342 11:27:19.46 +42:28:24.79 24.06 1.7664 + 0.0009 
NS_376 11:27:16.91  --42:28:30.62 22.69 0.6840 + 0.0001 
NS_407 11:27:06.62 +42:28:37.80 22.49 1.8553 + 0.0009 
NS.411 11:27:05.17 = --42:28:37.75 24.28 1.2456 + 0.0007 
NS_477 11:27:04.23  +42:28:56.36 23.00 1.8524 + 0.0009 
NS_481 11:27:13.66  +42:28:57.12 24.35 1.3236 + 0.0002 
NS_505 11:27:18.95  +42:29:06.64 24.40 1.2888 + 0.0002 
NS_511 11:27:18.69  +42:29:08.42 24.13 3.3255 + 0.0006 
NS_548 11:27:13.21  --42:29:30.75 23.14 0.8161 + 0.0003 
NS_610 11:27:00.97 +42:26:48.25 22.77 0.6628 + 0.0004 
N$.969 (Arc la)  11:27:13.92  +42:28:35.43 24.90 2.2355 + 0.0003 
N$S.1115 11:27:03.09 — 4-42:29:07.03 23.61 2.0585 + 0.0007 
NS_1477 11:27:16.65 — --42:27:23.57 21.04 0.7203 + 0.0003 


NOTE—NS numbers refer to the MSA slit identifications assigned when the obser- 
vations were designed. Positions are object positions as measured on NIRCam 


images. 


“The photometry will be presented by Pierel et al. (2023, in prep.) 


b This spectroscopic redshift will be presented by Chen et al. (2023, in prep.) 


€ The spectroscopic redshift for this arc was first measured by Polletta et al. (2023). 


corrects the atmosphere for ground-layer distortions via 
multiple artificial stars that are projected by laser beams 
mounted on each of the two 8.4m apertures (Rabien 
et al. 2019). This imaging achieves a mean point-source 
FWHM of 029 but a limiting magnitude of only 24.07 
(Frye et al. 2019). And finally, Spitzer Space Tele- 
scope/Infrared Array Camera (IRAC) 3.6 and 4.5 um 
images were acquired in 2017 (Cy 13, PID 13024, PI: 
Yan). The HST, LBT LUCI+ARGOS K-band obser- 
vations, and the Spitzer observations provide baseline 
imaging prior to the SN event but otherwise are super- 


seded by the deeper and higher-resolution NIRCam im- 
ages. 

Spectroscopy of G165 was acquired using LBT/LUCI 
long-slit during Director’s Discretionary Time on 2023 
Apr 16 (PI: M. Nonino) to measure the redshift of the 
SN host galaxy Arc 2. The 240"-long slit was aligned 
along two of the three images of this triply-imaged 
galaxy: Arcs 2a and 2b. The data were acquired using 
the G200 grating with the zJspec filter to cover wave- 
lengths 0.95-1.37 um. The observations were carried out 
in binocular mode resulting into a total exposure time of 
6 hours. The data reduction and analysis were presented 
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Figure 4. NIRSpec spectra of lensed sources at z ~ 1.78. Wavelengths are in the observed frame, and the ordinate shows P 
in units of 1077? ergs ^! cm? Åt. The G140M spectrum is plotted in green and the G235M spectrum in blue. Detected lines 
are marked. The images to the right of each panel show the respective source with MSA slit positions overlaid and are oriented 
the same as in Figure 1. Colors follow the prescription in Trilogy (Coe et al. 2012) with red showing F444W and F356W, green 
showing F277W and F200W, and blue showing F150W and F090W. The microshutter depicted in blue is the one from which the 
spectrum was extracted. The spectra are presented in order of star formation activity with more quiescent sources with weaker 
Ha emission lines and stronger 4000 À and Balmer breaks at the top to emission-line sources with multiple nebular emission 
lines at the bottom. These six sources uncover a diverse set of galaxy properties all contained in this single high-redshift galaxy 
overdensity. The spectrum for Arc G165.2a appears brighter and redder than the one for Arc G165.2c owing to the slit being 
better centered on the source. 
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Figure 5. NIRSpec spectra of Arc 1a (=Arc NS_969, top) and a nearby object (Arc NS_46, bottom). Wavelengths are in the 


observed frame, and the spectra are in Fy in units of 107 !? ergs^! cm? À-1. The G140M spectrum is plotted in green and the 
G235M spectrum in blue. Detected lines are marked. The images to the right of each panel show the respective source with 
MSA slit positions overlaid, and are oriented the same as in Figure 1. Colors follow the prescription in Trilogy (Coe et al. 2012) 
with red showing F200W, green showing F150W, and blue showing F090W. The objects were observed in the same triplet 
of MSA slits with Arc la in the easternmost segment and Arc NS_46 in the middle segment, neither perfectly centered in the 
respective segments. The blue outline in each image shows the microshutter from which the spectrum to its left was extracted. 


by Polletta et al. (2023). The spectra for Arc 2a and 2b 
both show stellar continuum, the 4000 À and Balmer 
breaks, and the [O II]A3727 emission line. These give 
the SN host redshift z = 1.783 + 0.002 (Polletta et al. 
2023). 'The redshifts of a few lower redshift galaxies 
intercepted by the long-slit were also measured which 
do not match the cluster redshift or that of any known 
image system and will appear in a different paper. 

MMT/Binospec long-slit spectroscopy was obtained 
in Director’s Discretionary Time on 2023 Apr 18 (PI: 
Willner) to measure the redshift of the brightest cluster 
galaxy (BCG) in the entire field, which resides in the 
northeast (NE) cluster component. The slit position 
angle 118°5 was selected to intercept the BCG and an- 
other cluster member in the southwest component that 
lacked a spectroscopic redshift. The data were acquired 
with the G270 grating to provide wavelength coverage 
from 3800 to 9200A. The total science exposure time 
was 1800 s. The observing conditions were good and 
reasonably stable, with atmospheric seeing ranging from 
1717-1725. 


'The data were reduced using the observatory pipeline, 
which performed the usual calibrations (bias correction, 
flat-fielding, wavelength calibration and relative flux cor- 
rection), as well as the coaddition of the three separate 
600 s exposures and extraction of the 1-D spectra*. Al- 
though redshift-fitting software was available, we opted 
to measure the spectroscopic redshifts using our own 
software that contains a library of spectroscopic features 
and a reference sky spectrum extracted from the pre- 
sky-subtracted data. A redshift was measured for the 
BCG of z = 0.3368. This redshift and others intercepted 
by the long-slit appear in Table 3. 

In all, the spectroscopy includes the following new con- 
tributions: 30 redshifts equating to 28 nonredundant 
individual sources (NIRSpec),2 redshifts equating to 1 
individual source (LBT/LUCI), and 6 redshifts corre- 
sponding to 5 nonredundant sources (MMT/Binospec). 
The total sums up to 35 new redshifts. When combined 


^ https: //bitbucket.org/chil_sai/binospec 
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Figure 6. Limiting magnitudes (5c) within a 02 aper- 
ture based on the corrected RMS maps (§3.1). The dashed 
vertical lines mark the mode of each distribution, and val- 
ues are given in each panel’s label. The dotted line marks 
the shallower peak, corresponding to the non-overlapping re- 
gions of Epochs 2 and 3. The deeper peak corresponds to 
the non-overlapping regions of the Epoch 1 images, how- 
ever the dithering causes the edges of the Epoch 1 images 
to contribute to the lower peak due to the lower exposure 
time in those regions. The tail of the deeper peak towards 
still fainter magnitudes results from the combination of the 
three epochs of imaging. It is not present in the F115W and 
F410M filters, which were acquired only in Epoch 1. 


with redshifts from the literature (Frye et al. 2019; Pas- 
cale et al. 2022), the spectroscopic catalog contains 175 
redshifts. Of these, 34 fall within +4000kms~! of the 


cluster mean redshift z = 0.348, and we consider these 
to be confirmed cluster members. Interestingly, there is 
a lack of detection of a relative velocity difference be- 
tween the NE and SW cluster components, implying a 
configuration oriented preferentially in the plane of the 
sky (Frye et al. 2019; Pascale et al. 2022). 


3. IMAGING RESULTS 
3.1. NIRCam Photometry 


The extraction of multi-band photometry broadly fol- 
lows the approaches of Merlin et al. (2022) and Paris 
et al. (2023), which balance the need to make faint im- 
age detections while also limiting the introduction of 
spurious sources. Initial source detection was performed 
using SExtractor (Bertin & Arnouts 1996) in a two-step 
HOT+COLD process similar to Galametz et al. (2013). 
'The JWST NIRCam F200W image was assigned as the 
reference image, which corresponds to the diffraction 
limit of the telescope. For severely blended objects, a 
separate catalog using F090W for detection was intro- 
duced owing to the sharper PSF. 

'The object fluxes and uncertainties were measured in 
each filter using APHOT (Merlin et al. 2019) by assign- 
ing Kron-like elliptical apertures, isophotal apertures, 
and circular apertures of diameter 0/3. To compute 
realistic photometric uncertainties, 5000 point sources 
from WebbPSF (Perrin et al. 2015) of known fluxes were 
injected into blank regions of the images, and fluxes 
and uncertainties were estimated using 071 apertures 
using APHOT with the RMS maps associated with the 
images. The RMS maps were then rescaled such that 
the RMS of the measured flux distribution was consis- 
tent with the values expected from APHOT. Because the 
mosaic was drizzled from 3 separate pointings, this anal- 
ysis was performed separately on each overlapping and 
non-overlapping region between the 3 epochs of data. 

'The final photometry is PSF-matched using PSF mod- 
els generated from WebbPSF, where all filters are de- 
graded to the PSF of the FA44W image using pypher. 
The F444W image is chosen because it has the largest 
PSF. The PSF models provided satisfactory convolu- 
tion kernels and are found to ameliorate the trend for 
the F200W — F444W colors to be bluer than their true 
values for multi-exposure image mosaics simulated for 
the CEERS project (Bagley et al. 2022). We refer to 
(Pascale et al. 2022) for details regarding its implemen- 
tation. 

The histograms of the limiting magnitudes were com- 
puted separately and are shown for each filter in Fig- 
ure 6. The double-peaked distributions demonstrate 
the multi-epoch experimental setup, with the PEARLS 
(Epoch 1) being deeper, and so showing a peak at 
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Table 3. MMT/Binospec Spectroscopy 


ID R.A. Decl. Mobs” Zsp 


BCG 11:27:16.70 +42:28:38.75 17.69 0.3368 
s101  11:27:14.05 +42:28:21.12 18.81 0.3427 
s102  11:27:08.04 +42:27:44.95 21.52 0.2757 
s103  11:27:20.51 +42:29:02.01 22.14 0.4135 
s104  11:27:35.83 +42:30:33.25 b 0.2300 
s105  11:27:32.70 +42:30:14.07 b 0.3509 


NoTE—The object identifications (IDs) are arbitrary 
and used only for convenience in this paper. Positions 
are object positions as measured on NIRCam images. 


“The AB magnitude is measured in the F090W filter 
that most overlaps with the MMT/Binospec spectro- 
Scopic coverage. 


b Source is outside the field of view of the JWST 
PEARLS program. 


higher limiting magnitudes per pixel. The tail of the 
Epoch 1 peak towards fainter magnitude emerges from 
the combination of all three epochs and is not present 
in the F115W and F410M filters, which have coverage 
only in Epoch 1. The data are shallowest in F115W 
(miim,Ap = 28.34) and are deeper in all of the LW fil- 
ters relative to the SW filters. 


3.2. SN HÜpe 


'Three new point-source images were identified in the 
Epoch 1 imaging close in projected proximity to Arc 2, 
and with an source plane separation estimated from our 
lens model of 1.5-2 kpc. The geometrical arrangement 
of the point source relative to the arc flips parity on 
crossing the critical curve between images Arc 2a and 
Arc 2b and again between the Arc 2b and Arc 2c im- 
ages, as predicted by lensing theory if the point source 
is associated with Arc 2. Since the point source follows 
the same lensing geometry as Arc 2 it is not a spurious 
source or low-redshift interloper, and so is likely at the 
same redshift as Arc 2. The triply imaged point source 
was easily bright enough (mpisow,aB = 23.91 mag) to 
have been detected in the earlier HST imaging had the 
source been present in 2016 (Figure 7). The geometry 
and transient behavior immediately suggested a super- 
nova in the Arc 2 galaxy, and this was confirmed by later 
observations as detailed below. For this paper, the three 
images (Figure 7) are referred to as SN 2a, 2b, and 2c, 
and we refer to 84 for a description of our lens model. 

Figure 8 shows the photometry from the Epoch 1 
imaging, with the reddest filters (F356W, F444W) re- 
moved due to contamination from the host galaxy. The 


Figure 7. Residual image (left panel) after subtracting 
HST/WFC3 F160W (middle panel) from JWST NIRCam 
F150W (right panel, degraded to match the F160W PSF). 
This subtraction led to the serendipitous discovery of SN 
HOpe in all three images of Arc 2. Arc 2b is brightest, fol- 
lowed by 2c and then 2a. The image order based on the 
predicted time delay from the lens model is 2a, 2c, and then 
2b. The HST/WFC3 imaging was obtained on 2016 May 30, 
and the JWST NIRCam PEARLS imaging was obtained on 
2023 March 30. 


light curve was applied only to the Epoch 1 photom- 
etry. The best-fit is for a Type Ia SN, where the light 
curve template was drawn from Hsiao et al. (2007). This 
classification predicts a second peak, which for z = 1.78 
will appear at observed wavelengths > 1.8m. The lens 
model predicts that the SN 2a image arrived first, and 
this image is faint in Epoch 1, potentially intercepted 
already following the the second peak. The model pre- 
dicts SN 2c should be the second image followed by SN 
2b, which apparently was at or very near the first peak. 
The light curve points, although sparse, were addition- 
ally best fit by a Type Ia SN model with >90% proba- 
bility, using the photometric light curve classifier from 
Rodney et al. (2014). These initial light curve fits are 
only approximations, and so are not useful for estimat- 
ing the time delay. At the same time, they were suffi- 
ciently compelling to have Epochs 2 and 3 approved and 
executed in our disruptive DDT program. 

'The photometry and spectroscopy in Epochs 2 and 
3 confirmed the Type Ia SN designation. The imag- 
ing in three representative filters that covers all three 
epochs is shown in Figure 9. As expected based on the 
lensing predictions, SN 2a is brightest in Epoch 1 in 
all filters and fades thereafter. SN 2b is the last im- 
age to arrive. It is relatively bright in all three epochs 
because all three are near its light-curve peak. SN 2c 
is the intermediate image, which is seen after the first 
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peak even during Epoch 1. Accordingly this image fades 
in bluer filters but remains bright especially in the LW 
filters. The full photometry, the methods for perform- 
ing the photometry corrected for background galaxy 
halo light and microlensing effects, and the photometric 
time-delay measurements, are presented by Pierel et al. 
(2023, in prep.). The supernova spectra, the SN type 
classification, and the spectroscropic time-delay mea- 
surement are presented in Chen et al. (2023, in prep.). 
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Figure 8. Photometry for SN HOpe based on Epoch 1 
imaging. The abscissas are in observed frame (lower) and 
rest-frame (upper), where the observation epoch was MJD = 
60033. The three images gave photometry for the three dif- 
ferent times sampled, with SN 2b arriving first (correspond- 
ing to MJD = 60033), followed by SN 2c, and then SN 2a, as 
indicated in the legend. Each point is further corrected for 
lensing magnification and time delays according to the best- 
fit lens model (§4). Points are slightly offset horizontally for 
clarity (such as for the F277W filter), and lines show best-fit 
model light curves color-coded by filter as indicated. 


3.3. Photometric Redshift Estimates 


Estimates of the photometric redshifts z,, were made 
across the grand mosaic using EAZY (Brammer et al. 
2008; Brammer 2021) and LePhare (Arnouts & Ibert 
2011). EAZY spectral energy distribution (SED) tem- 
plates were optimized for the identification of high- 
redshift galaxies in JWST/NIRCam imaging (Larson 
et al. 2022). A comparison of the 63 galaxies which have 
both spectroscopic and photometric redshifts shows 
good agreement for both approaches. In 95 % of cases, 
the photometric redshifts are within 15 % of the spec- 
troscopic redshifts. We briefly examine the only three 
outliers present across the two codes (Figure 10). The 
one source with a photometric redshift that is too high 
Zph ^v 6.3 — 6.4 compared to the spectroscopic redshift 
Zsp = 1.13 has a secondary redshift solution at zph = 1.5 
with a nearly equivalent goodness of fit. 


Of the two outliers with lower than expected photo- 
metric redshifts, one is found by LePhare at zph = 0.5 
with zp = 3.95, and one found by EAZY at zph = 0.08 
with Zsp = 3.3. The LePhare outlier is Arc 5.1a, which 
has a less favored secondary solution at zph ~ 4.1, a de- 
generacy which may be common at this redshift (Frye 
et al. 2023). The EAZY outlier has filter fluxes signifi- 
cantly enhanced by emission lines, LePhare's treatment 
of which is typically found to perform better, on average 
(Adams et al. 2023; Frye et al. 2023). 

Following this test of our photometric redshift ap- 
proach, we then extended the photometric redshift es- 
timates to the full multi-band object catalog. A pho- 
tometric redshift is considered secure if the object is: 
(1) in the field-of-view for all filters, (2) detected in a 
minimum of six filters, and (3) spatially resolved from 
its neighbors. The resulting distribution of photometric 
redshifts peaks at the cluster redshift and displays minor 
peaks at the z ~ 1.7 and z = 2.3 bins (Figure 10). The 
z = 1.7 bin corresponds to the redshift of the SN, pro- 
viding evidence that the SN host galaxy (Arc 2) may be 
one member of a larger galaxy group. Meanwhile, the 
peak at z — 2.3 aligns with the redshift of the lensed 
DSFG Arc 1, from which we infer that Arc 1 may be 
part of another background galaxy group. The physical 
properties of these two groups are discussed in 86.2. 


4. STRONG LENSING MODEL 
4.1. Lens Model Constraints 


The light-traces-mass (LTM) model requires as inputs 
the positions and masses of cluster members, the clus- 
ter redshift, and the image systems' identities, positions, 
and redshifts if known. We selected cluster members by 
their spectroscopic redshifts, when available, and aug- 
mented this list with sources selected by their near- 
infrared colors. In particular, the “1.6 um bump" is a 
feature of the stellar populations of >1 Gyr galaxies, ap- 
pearing when massive stars no longer dominate the com- 
posite galaxy emission spectrum (Sawicki 2002). For a 
cluster at z = 0.35, the 1.6 um bump appears as a posi- 
tive slope in the F090W — F150W color and a negative 
slope in the F277W — F444W color. Figure 11 shows the 
color-color selection for G165: F090W — F150W 0.5 
and F277W — FA44W « —0.5. The galaxies with spec- 
troscopic redshifts reassuringly occupy the expected re- 
gion of this color-color space. The cluster list is also in 
agreement with the “red cluster sequence" method for 
identifying cluster members (Gladders & Yee 2000), and 
enables longer-wavelength selection. The one outlier at 
F090W — F150W = 0.35 and F277W — F444W = —0.43 
has z = 0.3548 that places it in the cluster, but this 
galaxy’s photometry was skewed by nebular emission 
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Figure 9. Time series of SN HOpe imaging observations. Negative images in the three epochs are shown for each of the three 
images of the SN. Of the eight NIRCam filters acquired in Epoch 1 and six filters in each of Epochs 2 and 3, only three filters 
are depicted for clarity. The image stamps have the same orientation as in Figure 1 and are 3” on a side. The SN appeared 
first in image 2a, where it was intercepted on the decline after the second IR, peak in the light curve. SN 2a and SN 2b were 
intercepted near the first peak, and Epochs 2 and 3 also trace the second peak detected in the NIRCam LW channels of this 


Type Ia SN. 


lines. The main catalog of cluster members contains 
161 galaxies enclosed in a region of 3’ on a side centered 
on (R.A., Decl.) = (11:27:13.9143,+42:28:28.427) and 
depicted by the large green square in Figure 1. 

The initial lensing constraints include the 161 cluster 
members and the 11 known image systems (Frye et al. 
2019). New spectroscopic redshifts from NIRSpec in- 
creased the value of Arcs 5, 8, and 9 as lensing con- 
straints. A new lens model was then constructed and 
used to identify new image systems. Table A1 gives the 
complete list of image systems. 


4.2. Light-traces-mass Model 


'The LTM modeling approach is fast and allows for the 
straightforward incorporation of lensing evidence and es- 
pecially of the image multiplicities (Broadhurst et al. 
2005; Zitrin et al. 2015). The model presented in this 
study is the first NIRCam-based lens model, and builds 
and expands upon this previous work. We refer to Frye 
et al. (2019) and Pascale et al. (2022) for details regard- 
ing the construction and implementation of the LTM 
model in the G165 field. The base model for this study 
included the 11 known image systems as lensing con- 
straints. It also included the 161 cluster members iden- 
tified above, which were each assigned a mass propor- 
tional to the galaxy's luminosity with the mass-to-light 
ratio being a free parameter determined by the model. 

Starting with this base model, the redshifts for sys- 
tems 1, 2, 5, 8, and 9 were fixed at their spectroscopic 
values while allowing the parameters of the remaining 
image systems free to be fit by the model. Additional 
image systems were then introduced gradually, each 


time making sure that the fit improved. This process 
identified 10 new image systems which we apply as lens- 
ing constraints and which are reported in Table A1. For 
image systems with multiple components, we separated 
spatially resolved individual clumps, as recorded in the 
table. This is because the relative locations of individual 
clumps within an image provided additional constraints 
on the lens model. We stopped including these minor 
constraints when the x-square did not indicate an im- 
provement in the overall fit. 

Figure 12 shows the model for which x? is minimized. 
'This model reproduces the angular positions of input 
lensed images to an RMS difference of 0765. We note 
that while this uncertainty in predicting counterimage 
positions seems large compared to the NIRCam pixel 
size, it is not due to NIRCam WCS errors. Instead is 
due to the inability of the lens model to precisely locate 
all the lensing mass at the right locations in the model. 
Yet, the stated RMS is generally accurate enough to 
identify all the plausible counter images, and hence to 
refine the lens model. Based on our lens model, the 
lensing mass within (projected) 600 kpc of the cluster's 
luminosity-weighted center is (2.64-0.30) x 101^ Mo with 
the uncertainties computed by following the approach of 
(Pascale et al. 2022). This value is consistent with the 
masses estimated by Pascale et al. (2022) and by Frye 
et al. (2019) to within the uncertainties. 


5. SPECTROSCOPIC ANALYSIS OF THE ARC 1 
AND ARC 2 GALAXY GROUPS 


5.1. Diagnostic Tools 
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Figure 10. Photometric vs. spectroscopic redshifts. Points 
depict EAZY and LePhare redshifts as indicated in the leg- 
end. The panel on the right gives the histogram of pho- 
tometric redshifts using LePhare, which peaks first at the 
G165 cluster redshift and then again at the redshift of the 
SN of 1.78 and at the redshift of the DSFG (indicated by the 
dashed lines). Bottom panel shows |óz|/(1 + z) directly. 


'The NIRSpec spectra and NIRCam photometry are be 
modeled simultaneously to provide galaxy physical prop- 
erties and star forming activity levels. Modeling is based 
on FAST--4- (Schreiber et al. 2018; Kriek et al. 2009).° 
This SED-fitting code incorporated the NIRCam pho- 
tometry and all of the NIRSpec spectroscopic informa- 
tion (the G140M, G235M, and PRISM data), thereby 
accounting automatically for any strong emission-line 
features, which may not be apparent in the photom- 
etry alone. Models are based on a single burst with 
a range of decay rates 7 and a star formation rate 
SFR(t) x exp(—t/7), where t is the time since the on- 
set of star formation. The models that provide good 
fits to the JWST data return information relevant to 
this study such as the underlying stellar absorption, the 
dust extinction, and the stellar mass. For galaxies with- 
out Balmer emission lines, the dust extinction Ay was 
extracted from the FAST+-+ SED fit, and E(B—V) was 
computed based on the Calzetti et al. (2000) reddening 
law. 

Emission lines provide diagnostics on the galaxy clas- 
sification and activity levels (Baldwin et al. 1981; Curti 
et al. 2020, 2023; Li et al. 2023; Mingozzi et al. 


5 https://github.com/cschreib/fastpp 


2023, and references therein). The accessible emission 
lines are [O 1I] 43727, H8, [O II] AM4959,5007, Ha, 
[N II] AA6548,6584, and [S II] AA6717,6731. We used 
SPECUTILS to measure the line fluxes of Ha, H8, and 
other lines directly from the continuum-subtracted spec- 
trum or via multiple-line fits when Ha is blended with 
[N 11] A6584. These two line fluxes were corrected for 
underlying stellar absorption based on FAST++ esti- 
mates. Typical values for these equivalent width correc- 
tions were ~2 A for Ha and «4 Å for H£, both being in 
the expected ranges (Reddy et al. 2018; Sanders et al. 
2021). The Ha line flux was then corrected for dust ex- 
tinction based on the measured Balmer decrement (e.g., 
Domínguez et al. 2013). For Ha and Hó and for the 
computation of D(4000) we followed the Balogh et al. 
(1999) definitions. Otherwise, the line fluxes and equiva- 
lent widths typically were measured over a ~1000 km s^! 
width centered on the line core. Given the grating dis- 
persion of 6.4 À pix-!, this amounts to ~8 pixels. De- 
creasing this velocity width often resulted in a loss of sig- 
nal near the line tails. This study follows the usual con- 
vention that negative equivalent widths indicate emis- 
sion, and positive values indicate absorption, but ex- 
ceptions are that Ha and [O II] A3727 in emission are 
expressed as positive numbers. 

For galaxies that show strong emission lines, the Ha 
line flux and Balmer decrement give SFR and dust ex- 
tinction (Ay) directly. Other emission lines constrain 
the galaxy activity levels. The emission-line ratios 
O WI|/H6 (Ra), [N n]/Ha (N2), [S I]/Ha (S2), and 
O I]/Ha place the galaxy with respect to the Kew- 
ley et al. (2006) maximum-starburst line that separates 
star-forming galaxies (SFGs) from active galactic nuclei 
AGNs) (e.g., Mingozzi et al. 2023). The above line 
ratios plus ([O 111]|5007/H8)/([N 11]6584/Ha) (O3N2), 
[S 116716,6730/Ha)--([O 11]5007/H8) (R532), and 
[O in|5007/H8)/([S 1]6716,6730/Ha) (O3S2) allow 
multiple estimates of the gas-phase metallicity. 

In some cases the MSA slit does not fully cover the 
source galaxy's light, making the measured line fluxes 
lower limits. To compute a correction factor for incom- 
plete coverage of the lensed sources, we performed pho- 
tometry using the PRISM spectra, which provided un- 
interrupted coverage of all 8 NIRCam bands. This syn- 
thetic photometry is then compared with the NIRCam 
photometry, with the errors coming from the standard 
deviation of the differences between the 8 synthetic pho- 
tometric values from their true photometric values. This 
correction factor was then applied to the values in Ta- 
ble 4, where stated. On applying the corrections for un- 
derlying stellar absorption, dust extintction, and incom- 
plete slit coverage of the source, the Ho flux was used to 
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derive the SFR using the Kennicutt (1998) relation. The 
specific star formation rate (sSFR) was then computed 
by dividing the measured SFR by the stellar mass from 
FAST++. We note that this correction for incomplete 
slit coverage of the source are applied to the SFR but 
not to the stellar mass. This is because the stellar mass 
is computed by the SED fit scaled to the photometry 
integrated over the entire source. The uncertainties on 
the line fluxes and SFRs stem from the stated uncer- 
tainties of the flux values propagated by bootstrapping 
with a minimum of 100 realizations with the standard 
deviation of the distribution yielding the uncertainties. 
We proceed to investigate the physical properties of the 
Arc 1 and Arc 2 systems by a spectroscopic analysis. A 
more detailed spectroscopic line analysis that includes 
also radio information will appear in an upcoming pa- 
per (Foo et al. 2023, in prep.). 
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Figure 11. NIRCam color-color plot bracketing the rest- 
frame 1.6 um bump used to select G165 members at z = 
0.35. Small points show galaxy colors measured by NIR- 
Cam. Blue triangles indicate cluster members selected by 
spectroscopic redshifts, and green stars indicate Milky Way 
stars. The yellow-filled region in the lower right shows the 
color-selection area for cluster members. The colored lines 
shows colors of stellar populations of 3, 4, and 7 Gyrs at 
0.1 < z < 1.5 with redshift encoded as shown by the color 
bar. These colors represent a low-dust elliptical galaxy with 
no AGN. 


5.2. The Arc 1 System: Arcs 1a and NS 46 


The spectroscopically-confirmed members of the “Arc 
1 system” consist of two galaxies at z = 2.24. The spec- 
tra of Arc 1a (the DSFG) and Arc NS_46 present many of 
the emission lines expected of starbursting galaxies. For 


each object the ongoing SFR is measured off of the Ha 
line flux corrected for underlying stellar absorption and 
for dust extinction by the Balmer decrement method. 
The values corrected for lensing magnification are SFR= 
20173 Mo yr-! and 80729 Mo yr-! for Arcs la and Arc 
NS_46, respectively, where our lens model estimates u ~ 
2.7—2.9 for each. The SFRs are both high, and similar 
to within the uncertainties. Moreover, the two galaxies 
are closely-separated in radial velocity (71400 km s^ 1), 
making it tempting to ask if a galaxy interaction may 
explain the ongoing star formation. An estimate of the 
dust-obscured SFR, compared to unobscured (e.g. via 
Ha), may be ~50-90% (Whitaker et al. 2017). From 
which, we adopt a value of ~75% to represent the ob- 
scured, far-IR derived SFR to compare with our Ha- 
derived SFR. We suspect that the source does not have 
a 100% dust covering fraction because rest-frame UV 
light (observed-frame g-band) is detected (Frye et al. 
2019). This correction for obscured star formation leads 
to a rough estimate of the total SFR in Arcs la and 
NS_46 that is still higher by a factor of four. 


For Arc 1a, which is the lensed DSFG, a value for the 
SFR was measured encompassing the different counter- 
image (Arcs 1b/c), resulting in a total (rest-frame far- 
infrared luminosity) pS F R = 12,000—23,816 Mọ yr! 
(Harrington et al. 2016). The Ha vs far-IR-derived SFR 
values differ by up to ~order of magnitude when cor- 
rected for the lensing magnification factors obtained by 
our LTM model of ~2.9 and 30 for Arcs la and 1b/c, 
respectively. An attribute that makes their compari- 
son a challenge is that according to the lens model, Arc 
la corresponds to the total integrated starlight of the 
galaxy whereas Arcs 1b/1c equates only to a small por- 
tion of the galaxy light. Thus different regions of the Arc 
1 galaxy appear in Arc la and Arc 1b/c, respectively. 
Moreover, the light from Arc 1a is not easily separable 
from Arc NS.46 (Figure 5). Any cross-contamination 
of the light would alter the SFR tallies. A detailed 
lens model that separates out the independent contri- 
butions of Arc 1a from Arc NS.46 is beyond the scope 
of this study, and will be presented in a different paper 
(Kamieneski et al. 2023, in prep.). 

Arcs 1a and NS.46 are classified as star-forming galax- 
ies (SFGs) on account of their high Ha equivalent 
widths, > 10A (Li et al. 2023). Balmer family lines 
are detected in emission all the way up to and includ- 
ing Hd, and [O 11|A3727 is strongly-detected in emis- 
sion (Worry >5 A), thereby further classifying these 
two sources as short starbursting galaxies (Balogh et al. 
1999). On a larger timescale, the Hó absorption line 
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Table 4. Galaxies in Compact Groups at z ~ 2 
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ID R. A. Decl. Zsp Mpo0w E(B-V) u loguM* sSFR 
(J2000) (J2000) (AB) (Mo) (Gyr7!) 

Arc 2a 11:27:15.34  +42:28:41.05 1.783340.0010 — 20.30 0.22+40.04¢ 5.1 11.7  0.076+9-0407 
Arc 2c 11:27:15.98 +42:28:28.72  1.783440.0005 20.26 0.2040.04° 7.2 11.7 0.07079 0334 
NS.337  11:27:16.28 +42:28:23.59 1.7810 +0.0009 — 21.77 0.13+0.11° 2.5 10.7 gh 
NS.342  11:27:19.46  --42:28:24.79 1.7664 +0.0009 24.06 0.01+0.02° 16 9.18 2000 ee 
Arc 8.2c 11:27:15.80 +42:28:28.90 1.7839+0.0002 — 26.13 0.24+0.17° 9.0 8.45 Rotes 
NS.123 11:27:04.64 +42:27:15.24 1.7874 + 0.0003 23.94 — 0.32--0.05b 1.5 9.05 dst 
Arc9c — 11:27:16.11  --42:28:27.99 1.7816 +0.0002 — 26.10 — 0.312-0.05 ^ 5.0 — 7.43 2604120 
Arc la 11:27:13.92 +42:28:35.43 2.2355 + 0.0003 24.90 0.85+0.255 2.9 10.4 Eae 
NS.46 11:27:13.87  --42:28:35.67 2.2401+0.0002 24.04 0.81+0.04 2.7 9.51 orti 


NoTE— Column 1: lensed source name; Column 2: R. A.; Column 3: Decl.; Column 4: 


NIRSpec spectroscopic redshift; 


Column 5: F200W apparent magnitude; Column 6: Color excess due to reddening; Column 7: lensing magnification 
factor estimated from our lens model; Column 8: stellar mass as estimated from the FAST++ model, which is 
uncorrected for lensing magnification u; Column 9: specific star formation rate. 


? Masses are assigned a minimum uncertainty of 0.1 dex to account for systematics (Leja et al. 2019a,b). 


Measured directly from Balmer decrement 
© Derived from FAST++ SED fit 


d Corrected for incomplete slit coverage of the source. 


* Measurement of Ha line flux is not possible for this object. 


that typically appears in galaxies with star formation 
within the past 800 Myr is undetected in this data set 
(Goto 2007; Weibel et al. 2023). 

Taken together with its star-bursting galaxy classifi- 
cation, high E(B — V) value that is close to unity, and 
high metal enrichment levels, this galaxy pair overall 
shows some physical properties expected of a DSFG. 
For each galaxy in the pair the bulk of the star forma- 
tion is recent (< 5Myr), and possibly triggered by the 
interaction of Arc 1a with Arc NS.46. If this is the case 
then a galaxy-galaxy interaction may have also insti- 
gated any AGN activity. The measured and estimated 
internal properties of Arcs 1a and NS_46 are reported in 
Table 4. 


5.3. The Arc 2 system: Arc 2, NS.337, NS. 342, 
NS. 123, 8.2c, and 9c 


The *Arc 2 system" consists of Arc 2 (the SN host) 
and five other galaxies all at z — 1.78. Particular to 
Arcs 2a and 2c, nearly two dozen emission and absorp- 
tion line features are identified (Figure 4). Notably, for 
each image Ha is in emission while H£ is in absorption, 
making the Balmer decrement method of estimating 
dust extinction unavailable. Thus for these two galaxy 
images the Ho line flux was corrected for both underly- 
ing stellar absorption and for dust extinction from the 


FAST++ fit. The Ho line flux was further corrected 
for incomplete slit coverage of the source as described 
in 85.1. For Arcs 2a and 2c, we measured correction 
factors of 1.2 + 0.05 and 3.9 + 0.26, respectively. The 
SFR was then computed and then divided by the stellar 
mass. The stellar mass is (5.05 0.1) x 10!* Ma for both 
galaxy images, the highest of any other galaxy in this 
group at this redshift, while the sSFR has the lowest 
value of its cohort. 

The most conspicuous difference between the two 
galaxy images of this same lensed source is the high ap- 
parent brightness of Arc 2a (Figure 4). This is the case 
despite the fact that Arc 2c has a higher magnification 
factor. As can be seen in the MSA footprint (Figure 4, 
right-hand column) different regions of the galaxy are 
sampled in Arcs 2a and 2c, and in neither one does the 
slit cover the entire lensed source. In particular, the 
MSA slit for Arc 2a covers more of the central nuclear 
region of the arc, which may explain its higher bright- 
ness. Some clues may also be drawn from the rather dif- 
ferent values for D(4000) of 1.91 + 0.05 and 1.54 + 0.03 
for Arcs 2a and Arc 2c, respectively. It is interesting to 
ask if this 30 % difference may be a consequence of the 
MSA slit sampling a brighter central region that con- 
tains a higher fraction of stars in the emerging bulge 
galaxy component. 
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On galaxy classification, the Ha emission line indi- 
cates that at least some of the star formation is ongo- 
ing, and we measure Wiorr| > 5 A that corroborates this 
star-forming galaxy classification (Balogh et al. 1999). 
However, although Hó is detected in absorption, a fea- 
ture that is typically associated with galaxies which have 
undergone star formation within the past 800 Myr, it is 
weak (Wys < 5 À), possibly indicating that A- and early 
F-type stars are not yet dominating the spectrum (Goto 
2007; Weibel et al. 2023). There is a strong NaD line, 
which in combination with the higher D(4000) values 
indicate the presence of an older stellar population that 
formed within ~ 2Gyr (Wang & Lilly 2020). We note 
that the possibility of a post starburst (PSB) galaxy sub- 
category is generally ruled out by presence of Ha and 
[O 11]]A3727 in emission (Li et al. 2023, and references 
therein). 

The NaD absorption line is stronger in Arc 2 than 
for any other members in its z — 1.78 group, and may 
consist of both a stellar and an outflowing gas compo- 
nent (Cazzoli et al. 2016). Another prominent feature 
that is relatively rare at z ~ 2 is the detection of all 
three lines of CaT that is associated with cool stellar at- 
mospheres (Cenarro 2003). Of the culprits, supergiants 
would dominate the flux density budget, and may be 
another indicator of a previous star forming episode in 
this galaxy's more distant past (Z1 Gyr). In sum, Arcs 
2a and 2c are best-fit to a moderately dusty and mas- 
sive star-forming galaxy that appears to have already 
assembled much of its stellar mass by z — 1.78, a result 
that also matches the findings in Polletta et al. (2022). 

Arc NS_337 is apparently IR-bright (mr2o0w = 21.77 
AB mag). The prominent spectral features are detected 
in absorption except for He I, which is a common feature 
detected in nearly all of the z = 2 spectra in our sample. 
The computation of the SFR based on the Ha emission 
line is not available for this source. By its low W gs < 
5 À and Wiorr| <5 À, this galaxy is classified as passive, 
although the value of D(4000)—1.39 + 0.019 places it 
near the border with an SFG (Balogh et al. 1999). The 
addition of the low Wya < —3À along with the low 
Wiorn and Wps also disqualifies it as a PSB galaxy (Li 
et al. 2023, and references therein). 

Regarding its age, the Ha emission line is absent, 
implying at best only minimal ongoing star formation. 
Similar to Arc 2, Hó is detected in absorption, but is 
weak (Wys < 5 À), indicating that some star formation 
took place in the past ~800 Myr. There is a strong NaD 
line and a moderately-strong 4000 A continuum break 
D(4000) value, which indicate the presence of an older 
stellar population that formed within ~ 2 Gyr (Wang 
& Lilly 2020). Interestingly, similar to Arc 2 CaT is de- 


Figure 12. Color image using all NIRCam filters showing 
the central region of the G165 cluster. Magenta curves de- 
pict the z — 2 tangential critical curve obtained from our 
LTM model, and white labels show the image scale and ori- 
entation. Several giant arcs merge with the critical curve, 
an attribute that may assist with the detection of transients. 
'The lens model separates out the NE and SW cluster cores 
at this redshift. The color rendering is the same as Figure 3. 


tected, possibly uncovering an underlying population of 
giants/supergiants left over from an older epoch of star 
formation. The FAST4-4- fits estimate low dust extinc- 
tion that is 1 dex lower than Arc 2, and a stellar mass 
that is an order of magnitude lower than Arc 2. Overall, 
Arc NS_337 is a second relatively massive galaxy in this 
z — 1.78 group, and by a spectroscopic analysis appears 
to be quiescent. 

Arcs NS_342, NS_123, Arc 8c, and Arc 9c are emission- 
line galaxies with higher stellar activity levels, making 
them more akin to the z — 2.24 group spectra discussed 
in 85.2. Ha and Hf are detected in emission for most 
cases, enabling measurement of the SFR based off of 
the Ha line fluxes corrected for underlying stellar ab- 
sorption from a FAST+-+ fit, and for dust extinction by 
the Balmer decrement method. For Arc NS.342 where 
H£ is not detected in emission, we extract the value for 
Ay from our FAST--4- model fit. Relative to Arc 2, 
the sSFRs are all at least an order of magnitude higher, 
and the stellar masses are all ~ two orders of magnitude 
lower. 

'These four galaxies are all classified as SFGs based 
on the negative D(4000) values, and the high values of 
Wiorn > 5A and Wy_ > 10À (Balogh et al. 1999; 
Li et al. 2023). For Arc 9c the Balmer family lines 
are detected in emission all the way through Hô. For 
Arcs 8.2c and 9c we measure (Ws «0À) as well as 
(Worn 75 À), further classifying these two sources as 
short starbursting galaxies (Balogh et al. 1999). To get 
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their classifications by the optical diagrams, the flux ra- 
tios R3, No, S2 and [O 1]46300/Ho place Arc NS_342 
and Arc 9c in the star forming (SF) region, Arc 8.2c in 
the SF/composite region, and Arc N$.123 in the AGN 
(subclass Seyfert) region (Mingozzi et al. 2023). We note 
that only NS_123 shows evidence of having a predomi- 
nantly harder source of ionizing flux normally associated 
with AGNs. 

The presence of several rest-frame optical emission 
lines enables estimates of the gas-phase oxygen abun- 
dance. Ratios of emission line fluxes are computed for 
several of the line ratios and combinations of line ra- 
tios defined in 85.2 with the exception of Arc NS.342, 
whose diagnostic line ratios only include N2 and 5S5 since 
H8 is found to be in absorption even after accounting 
for underlying stellar absorption and correction dust ex- 
tinction. We obtain values for these four galaxies in the 
range of 12 + log (O/H) = 8.3—8.6. These sources have 
on average stellar masses and metallicities which place 
them higher than the mean relation correlating in Sloan 
Digital Sky Survey (SDSS) galaxies, the so-called “mean 
metallicity" relation (i.e., Curti et al. 2020). Taken to- 
gether with their high SFRs, it is interesting to ask if 
these counterparts of Arc 2 are intercepted during an 
epoch of a relatively-rapid build-up of stellar material. 
We refer to 86.2 for further discussion on the evolution- 
ary state of this compact galaxy group. 


6. DISCUSSION 
6.1. G165 Cluster Properties 


In a 1-D trace of the surface-mass density (Figure 13), 
the two peaks correspond to the cluster cores, and each 
core further separates out into more minor peaks that 
trace individual cluster members. An overall depres- 
sion appears between the two peaks in the 1-D trace, 
but not a sharp cutoff. We measure an angular sep- 
aration between the cluster cores of 38" equating to 
200 kpc at the redshift of the cluster. Inside of each 
core are two radio galaxies, with the two examples in 
the NE component further distinguished by showing ex- 
tended head-tail morphologies consisting of tails aligned 
~parallel to each other (Pascale et al. 2022, and ref- 
erences therein). Assuming a radial velocity of 600 km 
s^! typical of head-tail galaxies (Venkatesan et al. 1994), 
and that the radial velocity is similar in magnitude to 
the transverse velocity component, these quantities yield 
a crossing time of ~300 Myr. We note that the mean 
velocity difference of the two cluster core components is 
not well constrained given the current data set (Pascale 
et al. 2022). Redshifts of additional cluster members are 
needed to measure and/or better constrain the cluster 
velocity configuration. 


12° 24" 36" 48" 60" 72" 84" 96" 
Axis (arcseconds) 


Figure 13. G165 mass distribution from our lens model. 
Pink contours in the upper panel show the surface-mass (4) 
contours scaled to the critical value. The image is oriented as 
Figure 2. The lower panel depicts the 1-D mass distribution 
summed over 78" (—390 kpc) orthogonal to the collision axis. 
The value for & is summed up over the full field of view 
depicted in the diagram, which covers an angular range of 
96" (=480kpc). The NE and SW component mass peaks 
are marked with long-dashed vertical lines, and there is a 
modest depression between them. 


Pascale et al. (2022) report a relatively large off- 
set in the centroid position of the luminosity of 
the spectroscopically-confirmed cluster members which 
points to a major disturbance of the cluster. The ve- 
locity offset of the BCG from the systemic redshift of 
the cluster is another diagnostic of cluster virialization 
(Rumbaugh et al. 2018). The BCG’s measured red- 
shift of 0.3376 was originally drawn from the SDSS 
DR 17 archives and is independently-confirmed by the 
MMT/Binospec spectroscopy presented in this paper 
(Table 3). We herein confirm that both the brightest 
and the largest galaxy in the cluster entire (consider- 
ing both the NE and SW components) is significantly- 
blueshifted from the cluster's systemic velocity (~3400 
km s^!) that would place it near the outskirts of the 
cluster. 

One possible scenario is that the BCG is falling in to- 
wards the cluster from behind, such that the blueshift 
of the BCG and the redshift of the cluster are both 
~1700 km s-! relative to the velocity mean between 
the two of them. This relative newcomer to the cluster 
would not be detectable by our lens model as a subhalo 
within the NE component because the lens redshift is 
so similar. Moreover, the BCG is one of the two head- 
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tail radio sources, from which we infer that interactions 
with the interacluster medium are already imprinting 
the extended tails in the wake of its motion. This head- 
tail radio emission is less extended than the other clus- 
ter galaxy immediately to the north (see Pascale et al. 
(2022), their Figure 14). This possibly supports the view 
that the BCG is primarily falling along the line of sight 
towards the cluster rather than transversely, although it 
could also have a less-active nucleus. We note that for 
another merging cluster, El Gordo, the BCG was also 
significantly offset from the cluster’s systematic velocity 
( ~2400 km s^) (Frye et al. 2023). XMM/Newton ob- 
servations are approved for G165 (AO22, PID #92030, 
PI: Frye), which promise to offer some insights regard- 
ing its overall dynamical state by mapping out the mass 
distribution and searching for shocks. 


Figure 14. The Arc la galaxy group at z ~ 2.2. Photo- 
metric redshift fits identify a compact galaxy overdensity at 
z ~ 2.2 that surrounds the spectroscopically confirmed Arcs 
1a and NS_46 consisting of a total of 6 objects. These include 
the multiply imaged Arcs 3, 4, 6, and 15, as well a singly im- 
aged lensed galaxy located only 20 kpc away from Arc la and 
NS_46. Each galaxy exhibits a highly clumpy morphology 
consistent with cosmic-noon starbursts, and if their redshifts 
are correct, they may be interacting given their ~20 kpc 
physical separations (after demagnification according to our 
lens model). 


6.2. Compact Galaxy Groups & Associations 


Three peaks in the background galaxy distribution 
merit special attention (Figure 10). They correspond 
to: (1) a galaxy association centered in projection on 


the lensed DSFG Arc 1a at zsp = 2.24 (Figure 14), (2) a 
galaxy group centered at zp = 1.78 on the SN host Arc 
2 (Figure 3), and (3) a galaxy association centered at 
Zph = 1.65 located in the PEARLS G165 parallel field 
(Figure 15). These groups each have a different set of 
physical properties, as we discuss below. 

The Arc 1 Group: NIRSpec spectroscopy uncovers not 
just the DSFG (Arc 1a) but also a second galaxy proxi- 
mate to Arc la at a velocity separation of 1400 km s^! 
that appears to be the second member of an interacting 
galaxy pair (see Figure 5 and § 6.2). Such interactions 
are expected to be common at cosmic noon (e.g., Con- 
selice 2006), and may explain the star formation rates, 
corrected for lensing magnification of > 20 Mo yr^! for 
each one (Table 4). This galaxy pair is flanked by 5 other 
galaxies photometrically constrained to a similar red- 
shift of z ~ 2.2 (Figure 14). They include triply-imaged 
Arc systems 3, 4, 6, and 15, and a single galaxy im- 
age that we call LG1. These other galaxies have clumpy 
morphologies and blue colors consistent with active star- 
formation. Our lens model predicts for Arcs la, 3a, 4a, 
and 15a to be contained within a physical extent in the 
source plane of 20 kpc, and all of them to be contained 
within 40 kpc. 

If this tight configuration of galaxies is real, then com- 
bined with the evidence of starbursting activities in its 
two central members, it is possible that all seven of these 
galaxies may be engaged in interactions. The galaxies 
are expected to grow by star formation perhaps supplied 
by cold gas streams from the IGM (Dekel et al. 2009) 
and/or by galaxy mergers (Ellison et al. 2008; Scudder 
et al. 2012; Ellison et al. 2022). It is interesting to ask if 
we may be witnessing the “preprocessing” of the galaxy 
members, by which we mean that their stellar masses 
are built up prior to virialization and quenching (Ren- 
nehan et al. 2020; Sengupta et al. 2022). The answer 
to this question requires spectroscopy of the participat- 
ing members to extract their physical and star-forming 
properties, and so cannot be fully addressed with this 
data set. IFU spectroscopy in rest-frame UV and rest- 
frame optical would yield star formation rates and/or 
Lyman-a detections for all seven members that would 
confirm group membership and provide further insights 
into its evolutionary state. 

The Arc 2 Group: By far the dominant peak in the 
background galaxy redshift distribution is at z = 1.78 
(Fig 10). Hundreds of photometrically-selected galaxies 
with this redshift are distributed across the field of view 
of our NIRCam observations of ~1.5 x 2.5 Mpc at z = 
1.78 (Figure 10). This redshift is so common that it 
extends to about one-quarter of the image systems as 
well: Arcs 2, 8, 9, 10, and 16 (Table A1). If real, then 
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this galaxy overdense region effectively provides a grid of 
bright background z ~ 2 sources which may explain the 
abundance of giant arcs and image multiplicities in the 
field of this otherwise lower-mass example of a galaxy 
cluster. 

Of the six spectroscopically-confirmed galaxies at this 
redshift, our lens model estimates that the innermost 
four galaxies Arcs 2, 8, 9, and NS_337 labeled in Figure 3 
have a physical extent in the source plane of S33 kpc, 
with the SN host Arc 2 at the group center. These 
same four galaxies have a velocity spread of 900 km s71. 
The small spatial and radial velocity extents confirm 
the presence of a compact galaxy group at z = 1.78 
(Table 4). 

All but Arc 2 and Arc NS.337 have high sSFRs = 1 
Gyr-!. Arc 2 not only has a relatively low rate of ongo- 
ing star formation, but also shows little evidence of star 
formation in the past ~1Gyr (85.3). Moreover, Arc 2 
has a stellar mass that is 1—2 orders of magnitude higher 
than the other spectroscopically-confirmed members. It 
is interesting to ask if Arc 2 has already settled down 
after having built up the majority of its stellar mass, 
and is now surrounded by star-forming satellite dwarfs, 
potentially consistent with the picture of “downsizing” 
in hierarchical galaxy formation (Neistein et al. 2006; 
Fontanot et al. 2009; Oser et al. 2010). 

The z = 1.65 parallel field: The third group of galaxies 
is located at z ^ 1.65, and was uncovered in a photo- 
metric redshift search in the parallel field (Figure 15). 
These galaxies exhibit similar red colors, and on visual 
inspection all but one have apparent elliptical morpholo- 
gies. The morphology hints that the red color is a result 
of older stellar populations rather than dust, a scenario 
that is supported by our best-fit SED models. If so, then 
this galaxy association may consist mainly of quiescent 
galaxies. 

These galaxies have relatively high masses of 
1019711 M5 as estimated off of the SED models, stellar 
ages of Z1Gyr, and little to no dust extinction with the 
exception of the one SFG (yellow arrow in Figure 15). 
'This is consistent with the scenario of a build-up of stel- 
lar material at higher redshifts (z 2 3), followed by 
gas exhaustion at z S 1.5 just prior to an epoch char- 
acterized more by hierarchical growth (i.e., Rennehan 
et al. 2020). Spectroscopy is needed to confirm that this 
photometrically-identified galaxy association is a bona 
fide galaxy group, and may also yield better estimates 
of the SFR and dust extinction to corroborate this high- 
redshift quiescent galaxy group picture. 

Supernova Rate at z — 2: The total instantaneous 
(Ha-based) SFR in our NIRSpec sample that is cor- 
rected for lensing magnification (and dust obscuration 


Figure 15. The z = 1.65 galaxy group. Photometric red- 
shifts identified a group of seven galaxies within 15" + 20" 
in the G165 NIRCam parallel field. These galaxies broadly 
share similar colors and morphology, and photometric SED 
fitting finds all but one (yellow arrow) to be quiescent. These 
galaxies have a total physical extent of S130 kpc, making 
them another potential compact galaxy group. 


for Arcs la and NS_46) and integrated over both the 
Arc 1 and Arc 2 galaxy groups is 2500 Mo yr^!. AL 
though a rough correction has already been made to 
account for dust obscuration (85.2) it is likely to be con- 
servative, especially for the starbursting galaxies Arc la 
and Arc NS_46, which are especially dusty. In addition, 
there are other lensed galaxies at the same (photometric) 
redshift that could be star-forming but require spectro- 
scopic confirmation, such as Arcs 3, 4, 6, and 15. All of 
the above makes this SFR likely to be an underestimate 
of the true value. The SFR derived from the far-infrared 
luminosity Harrington et al. (2016) (corrected for mag- 
nification) of 400—800Mo yr~! uses spatially unresolved 
Herschel photometry covering Arc la and Arc 1b/c (in- 
cluding Arcs 3,4,6,15). At the same time only a portion 
of Arc la is imaged into Arc 1b/c. It should be empha- 
sized that a direct one-to-one comparison between these 
two SFRs derived by two different methods and applied 
to two different arcs is not entirely fair without near-IR 
spectroscopy of Arc 1b/c to serve as a calibration. 
Following the prescription in Young et al. (2008), and 
bracketing the two results (Ho and FIR-derived SFRs), 
we estimate a lower limit on the rate of core-collapse 
(CC) SNe 21-2 SNe yr^!. This value is consistent with 
results in ? in which CC SNe are summed up over the six 


JWST OBSERVATIONS OF G165 21 


Hubble Frontier Fields (Lotz et al. 2017). They found 
that with an JWST/NIRCam imaging cadence four 1- 
hour visits per year, the summed total rate integrated 
over all six HFFs is expected to be ~0.9 CC SNe and 
0.06 SNe Ia in F150W. Thus making regular monitor- 
ing of the G165 cluster alone may potentially be well- 
rewarded with the discovery of additional SNe. 


7. CONCLUSIONS AND FUTURE DIRECTIONS 


This overview study presented the full set of JWST 
observations in the field of G165. A transient was iden- 
tified in the PEARLS Epoch 1 NIRCam imaging relative 
to previous HST imaging (Figure 7). This point-source 
was situated 1.5-2kpc in the source plane from the mas- 
sive infrared-bright galaxy Arc 2 at z = 1.78 that is 
taken to be the SN host. Three different images of the 
Arc 2 plus SN pair were detected in the NIRCam imag- 
ing, and the pair flipped in parity with each crossing of 
the critical curve, as predicted by lensing theory. The 
three points had a fit to the light curve expected of a 
SN Type Ia (Figure 8). 

We identified 21 image systems which are used to con- 
struct the first NIRCam-based LTM lens model, which 
amounts to nearly a factor of two increase in the num- 
ber of image system constraints over previous works. 
Model inputs included 161 cluster members which were 
carefully selected by spectroscopy (where available), and 
the 1.6 um bump. On completing the first NIRCam- 
based lens model, we confirm the bimodal mass distri- 
bution. We also measure a mass within 600 kpc of (2.6 
+0.30) x 1014 Mo that is consistent with previous results 
also obtained by the LTM approach. 

NIRSpec spectroscopic redshifts were measured for 30 
lensed sources, which included all three appearances of 
the SN, two of the three images of Arc 2, image systems 
1, 5, 8, and 9, and other lensed sources. A spectrum of 
the BCG was also obtained using MMT/Binospec that 
confirmed it to be offset from the cluster systemic veloc- 
ity by ~3400 km s^ !, which we take to be evidence of a 
major cluster disturbance. The spectroscopic redshifts 
are presented in this paper (Tables 2 and 3). The focus 
of the spectroscopic analysis was on the galaxies at the 
same redshift as Arc 1 (z — 2.24) and Arc 2 (z — 1.78). 
These two systems of galaxies appear to be representa- 
tives of larger galaxy overdensities at these same red- 
shifts, an assertion that is backed up by minor peaks in 
the photometric redshift distributions. 

On the Arc 1 system, in a spectroscopic analysis and 
FAST++ fit to Arc la and its counterpart Arc NS_46, 
we measure relatively-high dust levels and high ongoing 
sSFRs expected of starbursting galaxies. By comput- 
ing the relevant line ratios across the rest-frame optical 


spectrum we find some evidence for AGN activity, and 
near-solar gas-phase oxygen abundances. Meanwhile, by 
other diagnostic absorption line features we infer the 
lack of significant star formation integrated over longer 
time scales, S1 Gyr. Our photometric redshift estimates 
pick out an additional five lensed sources characterized 
by clumpy star formation and situated within a spa- 
tial extent of 40 kpc in a group that includes Arcs 1a 
and NS.46. If real, then we speculate that this com- 
pact set of seven galaxies may be intercepted during an 
epoch of active star formation, a hypothesis that can 
be tested by obtaining additional and ideally IFU-based 
spectroscopy to characterize the star-formation in the 
various clumps/knots contained in each galaxy. 

Regarding the Arc 2 system, in a spectroscopic anal- 
ysis and a FAST++ fit to Arcs 2a, 2c and to its 
galaxy neighbors confirmed by NIRSpec which are Arcs 
NS_337, NS_342, 8.2c, 9c, and NS_123 (Figure 4), a 
stark contrast was unveiled between SN host Arc 2 and 
the other group members. Arc 2 dominates the stel- 
lar mass ((5.01 + 0.10) x 1014 Mọ), has a much lower 
sSFR (sSFRarc2a=0.0767 9-019 Gyr-!), and has only 
small amounts of ongoing star formation, which peters 
off to still lower amounts integrated out to ~1 Gyr. At 
the same time, Arcs 8.2c, 9c and NS.123 are emission- 
line sources which are on average 1—2 orders of mag- 
nitude less massive, yet have ~2 orders of magnitude 
higher values of ongoing SFRs scaled to their masses. 
Arcs 2, NS.337, 8.2c and 9c are separated in velocity 
7-900 km s-! and are situated within a spatial extent 
of 33 kpc. We infer that Arc 2 may be an example 
of a massive galaxy that has completed a star forma- 
tion episode ~1Gyr ago and that is now surrounded by 
star-forming satellite dwarfs, consistent with a galaxy 
downsizing scenario. 

The next steps focus on SN HOpe. The NIRSpec spec- 
troscopy of SN HOpe will be presented in an upcoming 
paper that confirms its classification as a Type Ia SN and 
measures a value for the spectroscopic time delay (Chen 
et al. 2023, in prep.). The full NIRCam photometry of 
the SN measured across all three observing epochs (and 
hence all 9 points on the light curve) that accounts for 
the subtraction of the Arc 2 galaxy halo light and cor- 
rection for any microlensing effects, and a measurement 
of the photometric value for the time delay will appear 
in a different paper (Pierel et al. 2023, in prep.). Mean- 
while, a more detailed lens model that performs a full 
surface brightness reconstruction of Arc 2 will take the 
next step in improving on the SN photometry (Caminha 
et al. 2023, in prep.). These efforts are complemented by 
a focused study on the dust extinction of Arc 2, which 
impacts understanding of the SN physics (Johansson et 
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al. 2023, prep.). A separate lens modeling investigation 
will focus on the compact galaxy group surrounding the 
DSFG Arc 1 (Kamieneski et al. 2023, in prep.), and an- 
other will focus on a more detailed spectroscopic analysis 
of the line features (Foo et al. 2023, in prep). A project 
is also underway to carry out a spectroscopic analysis 
of the quiescent lensed galaxy Arc NS274 (Valentino et 
al. 2023, in prep.). And finally, the time delay estimates 
generated from the lens models, the photometry, the 
spectroscopy, and the weighted combination of all time 
delay estimates from these models will be used to mea- 
sure a value for Ho (Pascale et al. 2023, in prep). Given 
the high rate of ongoing star formation across this clus- 
ter of 2400 Mo yr^1, regular monitoring of the G165 
field may be well-rewarded with the discovery of new 
SNe and other transients. 
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Table A1 lists the G165 arc systems used in the Section 4 lens model. The image-system designations for numbers 
1-11 follow Frye et al. (2019) with the exception of renaming Arc la to be the northernmost image and renaming 
Arc 1b/c to be the southern images. Ten image systems are new to this study. The table columns are: ID, Right 
Ascension, Declination, observed AB magnitude (SEXTRACTOR MAG_AUTO) in the F200W filter, spectroscopic redshift 
(zsp), photometric redshift (zpn), lens-model predicted redshift (zmoa), and the arc's discovery citation. Positions are 


J2000 on the GAIA DR3 system. 
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Table A1. Strong-Lensing Image Systems 


ID R.A. Decl. MF200W,obs Zsp Zph ^ Zmod Ref. 
la 11:27:13.91  +42:28:35.40 24.90 2.2355 + 0.00032 C18, F19 
1.1b  11:27:14.70 +42:28:23.81 22.18 2.2357 0d C15, H16 
1.16 11:27:14.80 +42:28:21.23 21.83 2.23570d C15, H16 
1.2b  11:27:14.74  --42:28:23.17 2.2357 0d C15, H16 
1.2c 11:27:14.78 — 4-42:28:22.12 2.23570d C15, H16 
2a 11:27:15.34 +42:28:41.05 20.30 1.7833 + 0.0010P€ C18, F19 
SNa  11:27:15.31 +42:28:41.02 G f g. This study 
2b 11:27:15.60 +42:28:34.19 19.89 1.783 + 0.002€ C18, F19 
SNb  11:27:15.60 +42:28:33.73 . f g. This study 
2c 11:27:15.988 +42:28:28.72 20.26 1.7834 + 0.00050 C18, F19 
SNc 11:27:15.94 +42:28:28.90 E Í g This study 
3.la  11:27:14.19 +42:28:32.15 2.21 C18, F19 
3.1b 11:27:14.24 +42:28:31.46 22.32 2.2 C18, F19 
3.1c 11:27:14.97 +42:28:17.37 C18, F19 
3.24  11:27:14.11  -F42:28:32.63 C18, F19 
3.2b  11:27:14.36 — --42:28:29.62 C18, F19 
4a 11:27:14.07 +42:28:32.66 22.96 2.25 2.173 C18, F19 
4b 11:27:14.37 +42:28:29.06 22.77 C18, F19 
4c 11:27:14.92 +42:28:17.29 23.62 C18, F19 
5.la  11:27:13.20 +42:28:25.73 25.36 3.9530 + 0.00042 F19 
5.1b  11:27:13.20  +42:28:24.62 24.88 F19 
5.1c 11:27:14.27 +42:28:08.83 26.14 F19 
6a 11:27:13.94 +42:28:32.75 24.62 2.18 2.215 C18, F19 
6b 11:27:14.36 — 4-42:28:27.64 24.32 C18, F19 
6c 11:27:14.84 +42:28:16.84 24.89 C18, F19 
Ta 11:27:15.24 +42:28:39.46 26.47 2.45 1.859 F19 
Tb 11:27:15.48 +42:28:33.71 24.13 F19 
Te 11:27:15.83 +42:28:28.33 24.23 F19 
8.la 11:27:15.27 +42:28:40.93 24.961 C18, F19 
8.1b  11:27:15.61 +42:28:33.04 22.99! C18, F19 
8.1c 11:27:15.87  +42:28:29.17 24.83^ C18, F19 
8.2a  11:27:15.28 +42:28:40.66 24.96 C18, F19 
8.2b  11:27:15.58 | 4-42:28:33.50 22.99h C18, F19 
8.2c 11:27:15.89  4-42:28:28.84 24.83! 1.7839 + 0.00025 C18, F19 
9a 11:27:15.43  4-42:28:40.57 25.54 C18, F19 
9b 11:27:15.57 . 4-42:28:36.02 25.33 C18 F19 
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Table A1 (continued) 


ID R.A. Decl. MF200W,obs Zsp Zph  2mod Ref. d 
9c 11:27:16.11 +42:28:28.01 26.10 1.7816 + 9.00022 C18, F19 
10a 1:27:15.18 -42:28:38.89 25.93 1.7 1.735 F19 
10b 11:27:15.46 -42:28:32.86 25.78 F19 
10c 11:27:15.74 +42:28:28.51 25.02 F19 
11.1a 1:27:15.77 -42:28:41.92 23.39 2.46 C18, F19 
11.1b 1:27:15.79 FA2:28:40.78 23.59 C18, F19 
11.1c  11:27:16.74 FA2:28:22.85 25.21^ C18, F19 
12a 1:27:14.52 42:28:23.48 27.61 2.581 This study 
12b 1:27:14.73 42:28:17.12 26.27 2.95 This study 
13.1a 1:27:14.04 42:28:37.76 26.86 4.209 This study 
13.2a  11:27:14.07 42:28:37.37 28.52 This study 
13.1b 1:27:14.75 42:28:27.59 27.89 This study 
13.2b 1:27:14.71 42:28:28.21 28.68 This study 
3.1c  11:27:15.14 42:28:18.82 28.89 This study 
13.20 11:27:15.16 42:28:18.38 28.80 This study 
4.la 1:27:13.09 42:28:26.83 24.62 2.25 2.342 This study 
14.1b  11:27:13.20 42:28:18.77 25.16 This study 
14.1c  11:27:13.84 42:28:11.32 24.69} This study 
4.2a 1:27:13.09 42:28:26.94 24.62? This study 
4.2b 1:27:13.21 42:28:18.51 24.21 This study 
14.2c  11:27:13.84 42:28:11.32 24.69 This study 
15a 1:27:14.11 42:28:34.10 26.38 2.339 This study 
15b  11:27:14.55 42:28:28.27 25.3 2.26 This study 
15c 11:27:14.96 42:28:19.09 26.75 This study 
16a 11:27:13.47 42:28:27.77 24.41 16 1.448 This study 
16b 1:27:13.97 42:28:23.90 25.20 This study 
16c 1:27:13.46 42:28:21.70 21.88 This study 
16d 11:27:14.22 42:28:13.49 24.51 This study 
17a 1:27:13.23 42:28:28.61 26.92 6 6.03 This study 
17b 1:27:13.22 42:28:23.00 26.89 This study 
17c 1:27:14.39 42:28:08.12 29.00 This study 
18a 11:27:15.35 42:28:37.70 26.27 This study 
18b 1:27:15.43 42:28:36.02 26.09 This study 
19a 11:27:14.18 42:28:39.44 27.26 This study 
19b 11:27:14.88 42:28:28.69 27.05 This study 
19c 11:27:15.30 42:28:20.11 27.10 This study 
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ID R.A. Decl. ™MF200W ,obs Zsp Zph Zmod Ref. d 

20a  11:27:14.44 . 42:28:35.63 - This study 
20b  11:27T:14.57 . 42:28:33.55 2 This study 
21a . 11:27:15.09 . 42:28:35.37 = This study 
21b  11:27:15.19 42:28:33.31 = This study 


?'lhe image system references are: C15 (Cahameras et al. 2015), H16 (Harrington et al. 2016), C18 


(Canameras et al. 2018), F19 (Frye et al. 2019) 


b Spectroscopic redshift for this arc was measured in this study. 


€ Spectroscopic redshift for this arc was measured by Cafiameras et al. (2015) or Harrington et al. (2016). 


d Spectroscopic redshift was measured at a different position along this arc. 


€ Spectroscopic redshift was first measured for this arc by Polletta et al. (2023). 


I The photometry will be presented elsewhere (Pierel et al. 2023, in prep.). 


I The spectroscopic redshift is assumed to be the same as that of the SN host galaxy Arc 2. The spec- 


troscopy will be presented by Chen et al. (2023, in prep.). 


A Photometry is blended with another image of the same system; the blended photometry is reported. 


